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ABSTRACT 


A study was conducted of the digital computer simulation 
of the equations of motion of the Surface Effect Ship with 
six degrees of freedom as developed by the Oceanics Corpora- 
tion. The version of the program under Sonenbeest ion was 2 
Simulation of the Bell Aerospace Systems 100-B Captured Air 
Bubble Craft (CAB) as adapted to the IBM-~360/67 digital 
computer and operating in irregular seas. The objective of 
the study was to optimize the simulation by reducing compu- 
tation time required to obtain solutions to the forces and 
moments acting on the craft while maintaining a reasonable 
degree of accuracy of the output variables. CPU time 
Savings achieved by use of the FORTRAN H compiler are 
documented. The dependence of CPU time and output accuracy 
on the tolerance levels established for the integration 
algorithm is discussed. CPU time savings versus output 


accuracy as the tolerance levels are increased is presented. 
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I. INTRODUCTION 


A. BACKGROUND 

In recent years the U. S. Navy has developed an interest 
in the Surface Effect Ship (SES) concept for combatant and 
Support type craft. As a result of this interest, Oceanics 
Incorporated was commissioned by the U. S. Navy to develop 
a digital computer simulation of the six degrees of freedom 
equations for SES craft that would yield time domain outputs 
of loads and motions. This initial report was delivered to 
the Navy Surface Effects Ships Project Office (SESPO) and 
dated August, 1971 [Ref. 1]. 

The Simulation program was designed to achieve a suffi- 
ciently accurate and stable soiution to the loads and motions 
equations for the modeling of the craft. The desire for real- 
time solutions was acknowledged, as they would be valuable 
for manned control Simulation studies and greatly reduce 
costly computer time during case studies [Ref. 1]. 

The program was developed using FORTRAN IV computer 
language and was prepared in modular form for ease in adapta- 
Llommmovartous SES configurations. The program delivered 
to the Naval Postgraduate School in October, 19/2, was a 
Simulation of the Bell Aerospace Systems 100-B Captured Air 
Bubble Craft (CAB), a 100 ton craft. 

Results of simulated runs conducted by Oceanics, 


utilizing a Control Data Corporation 6600 digital computer 


ec 


(using the scope 3.3 operating system) indicated that 
preplem Line to Computer time ratios, as low as 3:1 for 
calm water runs and as large as 1:4 for large amplitude, 
oblique, irregular sea cases at high speed, could be 
expected [Ref. 2]. 

The initial studies at the Naval Postgraduate School, 
using the IBM-360/67 computer with the standard Fortran G 
compiler to solve the 100-B loads and motion program for 
sea state conditions, resulted in problem time to computer 
time ratios of 1:40 to 1:70, depending on the type run 
simulated. As the sea state condition increased, the 
computer time increased to the point where for certain type 
fons, 2 to 3 hours of computer time were needed for problem 
solution [Ref. 3]. For example, broaching studies simulating 
65 knots with sea state 3 (irregular waves) at 150° enecunees 
angle and a 30 second run time, required 121 minutes of CPU 
time. It should be noted however, that in the broaching 
study, the output requirements included a listing of sidewall 
gap each time subroutine RHS is entered thus adding a signif- 


femme amount of CPU time. 


Bee Obs ECTIVES 

The high ratio of problem time to computer time for sea 
state operations demanded that an in depth study be made of 
the Oceanics SES simulation program as applied to the IBM- 
360/67 located at the Naval Postgraduate School's W.R. Church 


Computer Facility with the following objectives in mind: 


Ih: 





1.) Careful analysis of each subroutine to determine 
those areas where the bulk of the computer time is being 
spent. 

2.) Determine methods whereby the computational effi- 
clency of the program, can be maximized with minimum loss 
in accuracy. 

This thesis examines the SES 100-B digital computer 
simulation program on the IBM-360/67 computer with the 
above mentioned objectives in mind and makes recommendations 
for substantially reducing computer time while maintaining 
a reasonable degree of accuracy with regard to the loads 


and motion simulation. 
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II, INTEGRATION METHOD 


A. INTRODUCTION 

The SES simulation program is based on a mathematical 
description of the craft in six degrees of freedom as 
reported by Kaplan, Bentson, and Sargent [Ref. 1]. The 
resulting equations, in a form Suitable for numerical 


integration on a digital computer, are given by 


: F 
eee | " 
e F ; 
v= 2 - Ru ce) 
. +#F 
We — (3) 
pe kiten Oo Sn'txz (4) 
it Se ° «6t CUT I ¢ 
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Mm = Pa (Qs - one? Clee 


Where X, Y, and Z, are the 3-dimensional spatial coordinates; 
P, Q, and R, are the angular velocities about the X, Y, and 
Z, axis respectively; I is the mass moment of inertia, and 
F is the force. My is the mass of bubble air in the plenum; 
24 


bubble; Q 


- is the volumetric quantity of air flow rate into the 


nove is the flow rate due to leakage, and 9, is the 
Standard atmospheric reference density. 

The program is of modular design utilizing a parti- 
tioning technique on the craft to calculate the various 
forces and moments exerted on the hull as it moves through 
the water. The forces and moments acting on the craft are 
computed in various subroutines using initial conditions 
supplied by subroutine INCON. Subroutine RHS consists of 
a Simple summing algorithm to compute the total of the forces 
and moments in six degrees of freedom. 

Subroutine INTGRL calls subroutine RHS for the calcula- 
tion of the right hand side of equations defined by the 
Runge-Kutta—-Merson (RKM) algorithm for each required step 
size. The elements of each k vector in the RKM algorithm 
are the ten variables defined by equations (1) through (9) 


and equation (13). 
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piance the variables X, Y, and ) defined in equations 
(10), (11), and (12), are expected to have slow variation, 
it was decided that a simple means of integration using 
Simpson's Rule could be applied to these equations. These 
variables are calculated in the main program [Refs. 1 and 2]. 

A complete description of the computer program including 
its organization, flow charts, input and output description, 
program listing, and a user's manual for the operation of 


the program are available in Refs. 2 and 4. 


B. SUBROUTINE INTGRL 

The numerical integration technique used in subroutine 
INTGRL is the Runge-Kutta-Merson (RKM) method for solution 
of a system of eee order differential equations of the 


form 
y = f(t,y) (14) 


the recursion formula is 


k. t4Uk, tk 
- il, y 5 8) 
es ee 5) + o(h~) , C25) 
where 
ae: 
dei h 
Ee = 3 f(t, +3 ; Yi tk) CL) 
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ky k 


h 2 


Se 
Ke " z F(t, +s , %, t isle — (18) 
k 9k 
oh ok 3 
Ky peo aie, ee ge he) (19) 
3k Ok 
eel bi a5 
ke = 3 f(tpth , Y, + = - 42 + 6k, ) (20) 


and h is the stepsize. The truncation error in (15) is 
given by 


9k k 
i 3 5 
E—E = kK = + AKy — oe (21) 


1 2 

Initial conditions are read in through subroutine INCON 
to establish the starting point for each of the ten variables 
plus a trial stepsize (h) and a maximum error tolerance level 
for each integrator. Computations are then carried out and 
the truncation error given by equation (21) is computed. If 
any element in the error vector € is greater than the corre- 
Sponding maximum tolerance level then the trial time stepwise 
is halved and the computations are repeated, If the error is 
less than the tolerance level, but greater than 1/16 this 
level, the computations proceed in accordance with the RKM 
algorithm. If all elements in e€ are less than or equal to 
1/16 the maximum Peete orice level, the stepsize is doubled 
in the next integration step. This procedure automatically 
adjusts the time step in the computation to reflect the 


nature of the perturbations being experienced by the craft. 
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As a result, for slowly-varying phenomena (e.g., calm 

water, slow speed simulation), the computations are carried 
out quickly since larger time steps are taken, while high 
frequency effects (e.g., heavy seas, high speed simulations), 


will cause an increase in computation time. 
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Lit. ‘PROCEDURE 


A. INTRODUCTION 

The obvious starting point for an efficiency study of 
the SES simulation program was to gain a thorough knowledge 
of the force and moment algorithm used. The object was to 
Study the interconnections and interactions between the 
various subroutines with a goal of arriving at a more effi- 
cient (i.e. less time consuming) solution by reprogramming 
and/or better utilization of the IBM-360 software. 

A major reprogramming effort was considered to be beyond 
the oe of this thesis. However, investigation of subrou- 
tines which were known to consume the bulk of the computation 
time was conducted. In conjunction with this investigation 
various timer functions built into the IBM-360 were utilized 


CO pin Pelintetime consuming procedures. 


B. IBM-360 COMPILER 

The SES simulation program was originally adapted to the 
Naval Postgraduate School IBM-360 computer utilizing the 
Standard FORTRAN G compiler. This compiler features a DEBUG 
facility and core optimizations which lends itself to the 
Smaller job, teaching environment prevalent at the school 
(Ref. 5]. Among other compilers, the IBM-360 computer system 
has available FORTRAN H, an IBM processor which produces a 


huehaly optimized code which is highly desirable for large 
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production programs such as the SES simulation program 
(Ref. 6]. 

Compilation of the program in FORTRAN H resulted in an 
average CPU time savings of thirty-seven percent over the 
FORTRAN G compiler while yielding exactly the same output 
variable values. For example, a simulation of a thirty 
second run in a state three, head sea compiled in FORTRAN G 
required 80.5 minutes of CPU time. The same simulated run 
mine compiled in FORTRAN H required only 50.7 minutes of 
CrUstame. Avnumber of runSj under varying inivial ceondayaens- 
were made to verify the CPU time savings using the FORTRAN’ 
H compiler. Unless otherwise indicated all further runs 
discussed in this thesis were made using the FORTRAN H 


compiler. 


C. TIMER SUBROUTINES 

The IBM-360/67 at the Naval Postgraduate School's W.R. 
Church Computer Facility has several built in timer subrou- 
tines which were helpful in determining those areas in the 
program in which the bulk of the CPU time was being used. 
Two such routines were used in this study. 

1. Subroutine ILONUM 

Subroutine IONUM is a program which monitors the 

number of times an input or output device is accessed during 
the run of a computer program. In adapting the Oceanics 
simulation program to the IBM-~360, extensive use of discs 


as temporary storage devices was made. By monitoring the 


a 





number of input/output operations with IONUM a more effi- 
clent data blocksize for the dises could be determined, 
thus reducing the number of accesses required. 

es *ubrewtpines PROGHOOK 

several areas of the SES simulation program were 
suspected to be inefficient and/or requiring excessive time 
during the program solution. In aeciens to pin point these 
areas of the program, subroutine PROGLOOK was used. 

PROGLOOK, when linked to the program under consid- 
eration, monitors the sequential flow of the program instruc- 
tions. The routine "strobes" the instruction flow at a rate 
of 50 samples per second or approximately once every 4000 
instructions. In this manner a time history of the program 
functions was obtained. The output of subroutine PROGLOOK, 
which was of interest to this study, is in the form of a 
histogram which gives a plot of percentage of time the program 
spent performing various functions. Correlation of the loca- 
tion of these functions with a computer map listing of the 
SES simulation program served to point out time consuming 


areas within the program. 


D. INTEGRATOR TOLERANCES 

A prime suspect area of the simulation program for excess 
time consumption was the Runge-Kutta-Merson (REM) aleortenm 
used to solve 10 of the force and moment differential equa- 
fions. if during a given time increment, any one of the 10 


error function values is found to be outside the given 


ie 





tolerance, the stepsize is halved and all variables are 
computed using the new stepsize. This process is continued 
until all 10 error functions are within their respective 
tolerances and then a new time increment is started. 

The above stated procedure does not constitute an inord- 
inate amount of computer time under calm water simulations 
as the variable values tend to vary slowly and the RKM 
algorithm converges rapidly to the established tolerance 
values. However, craft perturbations caused by such condi- 
tions as simulations of rapid turns, irregular waves to 
simulate sea State, and 'oss of a propulsion Uityeeomecace 
some of the variables to vary at a high rate of change, thus 
requiring many iterations of the algorithm for convergence 
to the established tolerance levels and resulting in large 
CPU times. 

Reference 3 established the standard maximum tolerance 
levels for each integrator according to procedures set forth 
in Ref. 2. These tolerance levels were considered to be 
the basis of accuracy desired for the SES simulation program 
used at the Naval Postgraduate School and are listed in 
Table I. The input and output variables for ais of the 


13 integrators are shown in Fig.l. 
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Table I 


Tolerance Values For Integrators 


Integrator Number Output Variables Tolerance Values 
1 ~0001 
BUCO 
.000001 
.0000001 
~OOma 
.000001 


.000000001 


Do © yy © Be = “= Fe 


.000000001 


NS 


~0001 


io an n Ww LFS W WN 


~O001 


o 


1. Steady State Conditions 

‘A number of runs were conducted initially, for each 
simulation study made, in order to establish steady state 
conditions for a given set of initial conditions. Reasonably 
accurate steady state conditions are essential to prevent 
initial transients in the simulation program which can 
cause the program to abort, or worse yet, to give erroneous 
results. The key variables to initialize for a given input 
are thrust, draft, plenum pressure, and to a lesser degree 


Patem angle. 
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RKM  Runge-Kutta-Merson Method 
SR Simpson's Rule 


Bigure 1, Inputs and Outputs of Integrators 
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Two run conditions were selected to fulfill the 
objectives of this thesis. (1) Sea state 3, 40 knot speed, 
head seas, and (2) Sea state 3, 60 knot speed with turn. 
Table II gives the initial conditions for these two run 


eend2tions. 


Table II 
Initial Conditions 


(Sea State 3) 


Speed brat t Plenum Thrust ricen 
(kts) (In) Pressure (one Engine) Angle 
CES) (Lbs) (Deg) 

40.0 30.0 05. 10) TODO eo 0.50 


60.0 28.0 80.0 10800.0 0.45 


2. integrators 3 and 10 

The greatest effect of the sea state condition on 
the SES program computation occurs on integrator 3 (heave 
accelerations) and integrator 10 (time rate of change of 
the bubble air mass). The majority of the forces and 
moments acting on the craft when encountering sea states 
is manifested in changes in Z directed accelerations and 
plenum volume (i.e., bubble mass). 

hs a first approach to reducing the CPU time for 
sea state simulations, the changing of the implementation 


of the RKM algorithm was considered. It was felt that if 
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Integrator 3 and 10 could be isolated and their solutions 
made to converge before the other eight integrations were 
attempted, then a considerable computation time savings would 
be gained due to the reduced number of iterations needed 
through subroutines INTGRL and RHS. However,,due to the 
functional dependence of heave acceleration and bubble mass 
on some of the variables, it was found that Integrators 3 
and 10 could not be made to converge independently of the 
Other eight. 

Since both Integrator 3 and 10 were extremely sensi- 
tive to craft perturbations the decision was made to loosen 
the tolerances on these two integrators and compare the 
results with those obtained using the standard tolerances 
(See Table I). It was felt that if the percentage change 
in variable values was not too great (less than 10 percent), 
acceptable results could be obtained with a considerable 
CPU time savings. 

The two integrator tolerance levels were increased, 
in tandom, by a factor of ten in steps of one for a total 
of ten runs under a given set of initial conditions. Each 
computer run simulated a 30 second craft run and parameter 
values were printed out every 0.05 seconds for a sample 
rate of 20 per second per variable. In addition, the root 
mean square (RMS) value and the average value of the 


parameters was calculated over the entire run interval for 


ai 





comparison with values obtained from runs uSing the standard 
tolerances. Using these results an error analysis could 


be conducted. 


3. Integrators 1,2 and 4-9 

Analysis of the CPU time versus output accuracy 
as observed in the above runs indicated an optimum tolerance 
or 0.000007 Tor integrator seana 0. C0CZeror integrator 10. 
It was then decided to investigate the sensitivity of the 
other eight integrators. The method used for this investi- 
gation was to hold the above tolerances on integrators 3 
and 10 while increasing one of the other integrator toler- 
ances by a factor of ten (the other seven tolerances remain- 
ing at standard values). Each integrator was tested in 


Succession and the resultsS compared to standard runs. 
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IV. DISCUSSION AND EVALUATION OF RESULTS 


A. FORTRAN H COMPILER 

The efficiency of the FORTRAN H compiler can best be 
illustrated by the following example. Reference 3 cited 
the disparity in CPU times between simulation runs using 
control statement 00201 and 00202. These runs were made 
using the FORTRAN G compiler. Control statement 00201 inputs 
directly to subroutine INCON the overall weight of the craft 
plus the location of the center of gravity and the mass 
moment of inertia about the X, Y, Z, and XZ-axis. Control 
Statement 00202 inputs distributed weights and moments and 
Subroutine INCON must then compute the craft weight, CG 
and moments. It would appear that the additional computa- 
Wionemaivorlved Ieusing control! Statement O0202 would require 
mome CPU time than statement 00201. 

Results of Ref. 3 indicated that the opposite result 
was obtained and that CPU time decreased by approximately 
10 percent when control statement 00202 was used as input. 
Mie reason for this result has not been established and 
Should be looked into by future investigators. Nonetheless, 
results of runs compiled in FORTRAN H reduced the CPU time 
difference cited in Ref. 3 to less than 1 percent which 


Liimstrates the optimization capabilities of this compiler. 
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B. INTEGRATOR TOLERANCES 

The initial runs were made to verify the sensitivity 
of the Runge-Kutta-Merson algorithm to integrator tolerances. 
The runs simulated 40 knots speed, state 3 seas, head waves, 
SU }secona run, and confirmed that integrators’ 3 and 10 
were the most sensitive to craft perturbations. The results 


of these runs were as follows. 


Table III 
CPU Time Comparison 


(Sea State 3) 


Integrator Tolerance CPU Time PERCENTAGE Decrease 
..000001* 057,0.0,0.).> i835 5 Standard 
0.000002 0.0002 38.74 27.66 
0.000003 0.0003 Bo) oe, 3 Tied 
0.000004 0.0004 5136 41.40 
0.000005 OL TONE ONS 29.20 45.47 
0.000006 0.0006 28.90 4S . 73 
0.000007 0.0007 26.99 49.56 
0.000008 0.0008 26.19 51.09 
0.000009 0.0009 ap 03 Salo 
0.000010 0.0010 24.50 54.25 


* Standard for comparison 
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ie Integrator 3 and 10 


The tolerances on integrators 3 and 10 were loosened 
during the above stated runs while the other eight integra- 
tor tolerances were maintained at their standard values. 
Table Iii shows the results of these runs with regard to 
CPU time and Figure 2 is a graphical representation of the 
data. 

As the integrator tolerances were varied the RMS 
and average values of the output variables were compared 
with those obtained using standard integrator tolerances. 
The percentage deviations from standard values were computed 
in order to have a measure of output variable accuracy as 
a function of CPU time savings. Table IV tabulates these 
mesyits for maximum to@erance error levels. 

As was expected, the highest percentage deviations 
occured with heave acceleration (integrator 3) and the 
variables concerned with air flow (integrator 10). All 
other variables had insignificant deviations (less than 12) 
from standard values. In general, the percentage deviation 
of the output variables increased almost linearly as the 
integrator tolerances were loosened, with the maximum 
deviation occurring for maximum tolerance level. 

Increasing the Integrator tolerance levels above a 
factor of seven did not greatly decrease the CPU time required 
forepropvem solution but did tend to decrease the accuracy 


Peeriemerteeved OULDUL variables. For example, the initial 
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Output Variable Maximum Deviation 


From Standard Values 


Output Variable 


RMS 
Draft 0.04 
Pitch Angle 0.82 
Plenum Pressure Oo. 44 
CG Heave Acceleration eee) 
surge Speed G03 
X-Displacement 0.04 
Fan Power ab 
Air Flow Into Plenum O15 OFS) 
Air Flow Out of Plenum Caer 8 
Bubble Drag 0.61 
Pitch Rate O00 


tolerance increase of a factor of seven gives a CPU time 
savings of 26.56 minutes, while increasing the tolerance 


level to a factor of ten achieves only another 2.49 minutes 


Ole 


SS ££: oo oS 


Cy elo. a= 


Sn =] ai 


Maximum Deviations from 
Standard Value (Percentage) 


Average 


07 


. 00 
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00 
a¢2 
O04 
»33 
75 
nee 
By 
00 


savings in CPU time. At the same time the percentage 


deviation in the value of heave acceleration increases from 


Tish Seton am 
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2. Integrators 1,2 and 4-9 


With integrator 3 and 10 tolerances at 0.000007 and 
0.0007 respectively, each of the other integrator tolerances 
were increased in turn by a factor of ten. Simulated run 
conditions were as previously described and CPU times were 
compared. No appreciable CPU time savings was observed 
during any of these runs. The average CPU time was 27.59 
minutes for these runs. The percentage deviation in the 
values of the output variable was negligibly small, indi- 
cating that these variables were much less susceptible to 
the moments and forces caused by the simulated sea state 


Ponda >Lons. 


Ce DIFFERENT TYPE RUN CONDITIONS 

Various run conditions were simulated in order to study 
the effect of loosening the tolerances on integrator 3 and 
10 while maintaining standard tolerances on the other eight 
integrators. Each computer run Simulated a 30 second craft 
iene, 

Table V shows the results of these runs with respect to 
CPU time savings. 

Comparison of the resultS Shown in Table III and Table 
V indicate that the percentage of CPU time savings is about 
the same regardless of the type run simulated. Therefore 
increasing the tolerance levels on integrators 3 and 10 
by a given amount will result in a predictable amount of 


computer CPU time savings. 
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Table V 
CPU Time Comparison For 


Different Type Run Conditions 
(Sea State 3) 


Speed Integrator Tolerance Wave Encounter CPU Fereentare 
Asap: Decrease in 
(mots) No. 3 No. 10 Angle in.) CPU tire 
40 0.000001* 0.0001% ESS be.U9 Standard 
0.000007 0.0007 24,73 Serpe 
0; 000010 0.0010 23.47 54.94 
60 0.000001* 0.0001* Partial Turn 68.04 Standard 
0.000007 0.0007 34.34 4g,2u 
0.000010 0.0010 BIS SS Dono 
60 0.000001* 0.0001* 180° Turn 69.97 Standard 
0.000007 0.0007 36.36 48.03 
0.000010 0.0010 32.46 53.61 


* Standard for comparison 


In comparing results of deviations of output variables 
under different run conditions, no change was observed 
between the 40 knot, head sea simulation and the 40 knot 
150° wave encounter simulation outputs. However, under the 
more severe conditions of 60 knots with turns applied some 
of the output variable values did have an increased deviaticn 
from the standard values. Table VI depicts these results 


for the variables concerned. 
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Table VI 
Output Variable Maximum Deviation 


From Standard Values 
(Sea State 3, 60 knots) 


Output Variable AN, |G@pguelal Palen Maximum Deviation From 
Standard Value (Percentage) 
RMS Average 
Plenum Pressure Partial LAS 0.29 
7O0° 1330 p05 
CG Heave lesevigk Galati Mas nF 10\3) 
Acceleration 
Leg Spee: 0.00 
Fan Power Partial 0.90 4.58 
Poor 1.41 Mo. 
Air Flow Out Partial Oe 7 ineoo 


of Plenum 


D. SUBROUTINE PROGLOOK 

In order to determine other areas of the SES simulation 
program where large amounts of CPU time was being used, 
Subroutine PROGLOOK was utilized. Integrator 3 and 10 
tolerances were 0.000007 and 0.0007 respectively. A 40 
knot, state 3 sea, head waves, 30 second run was simulated. 
Piewve ssels @ histogram of the results of thé run. 

It was found that 41% of the computing time was spent 
imeeubrouvine WAVES. The bulk of this time was consumed 


in the wave table. A surprisingly large amount of time 
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(22%) was spent computing trigonometric function values. 
The bulk of the time spent in a given Subroutine was gen- 
erally found to be involved in repetitious DO LOOPs. 

The area labeled "Other" included various FORTRAN built 
in functions and the remaining subroutines and functions 
of the SES program which were not identified on the graph. 


Each of these routines contributed less than 1% of computation 


time. 
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V. CONCLUSION AND RECOMMENDATIONS 


A. INTEGRATOR TOLERANCES 

The computer run time dependence of the SES simulation 
program on the tolerance levels chosen for the Runge-Kutta- 
Merson algorithm is clearly demonstrated by the results 
shown in this thesis. Further, CPU time is directly 
related to the magnitude of the tolerance levels of inte- 
grators 3 and 10 and is relatively independent of the 
tolerance levels established for the other eight integrators. 

Output variable accuracy as would be expected is also 
related to the tolerance levels of the integrators, however 
the percentage deviations from standard values is small and 
the values obtained could be considered adequate Tor most 
purposes. In particular, initial studies under a given set 
of conditions would realize a considerable amount of CPU 
time saving by uSing increased tolerance levels, while 
maintaining sufficiently accurate solutions to determine 
the relative magnitude of the forces and moments acting on 
the craft. AS more precise values are required, the tolerance 
levels could be progressively tightened with an accompanying 
predictable increase in CPU time requirements (a valuable 
planning tool). 

It should be noted that the individual programmer should 
determine his standard of accuracy for each output variable. 


Since the magnitude of the outputs vary greatly, a small 


oe 





percentage deviation in one value may be acceptable for a 
given case study while being entirely too large for another 
output value under the same conditions. For example, the 
air flow rate into the plenum ‘has a magnitude on the order 
of 5000 cubic feet per second with a 1% deviation from 
standard which may be entirely acceptable, whereas the CG 
heave acceleration magnitude of 0.32 "g" units with 1% 
deviation may not be acceptable. 

Due to time limitations, the simulated runs conducted 
in this thesis were confined to relatively slow speeds in 
the 40-60 knot range. The deviations observed for the 
critical output variables of heave acceleration, plenum 
pressure, and air flow in and out of the plenum, increased 
with an increase in speed and with the introduction of turns. 
Therefore it is recommended that further studies be made at 
higher speeds and more severe sea state conditions to deter- 
mine the validity of results obtained with increased 


tolerance levels. 


B. DIGITAL PROGRAMMING TECHNIQUES 

Time studies of the SES Loads and Motion program uSing 
subroutine PROGLOOK indicated that an appreciable amount of 
computer time was spent in the trigonometric functions. Most 
large computers use nested Chebyshev Polynomials or Taylor 
Series expansions as an efficient and accurate solution to 
the trigonometric functions [Ref. 6]. A faster method of 


solution, for example, would be a four place table lookup 


4Q 





for the functions. A careful analysis of the loss of 
accuracy incurred by this method would have to be conducted. 
Another consideration would be the increased core require- 
ments due to the addition of a lengthy table. 

It was noted that considerable computation time is spent 
in the execution of DO LOOPs throughout the program. It is 
recommended that an in depth study be made with the goal of 
replacing the iterative DO LOOP process with straight forward 
calculations where possible. An excellent example of this 
type of programming is presented in Appendix A of this thesis. 
A straight forward programming technique resulted in the 


elimination not only of DO LOOPs but a Subroutine as well. 


C. HYBRID COMPUTATIONS 

Plans are being formulated at the Naval Postgraduate 
School with regards to the feasibility of adapting the SES 
Loads and motion simulation programs to a hybrid computation 
technique. The development of linearized equations of 
motion presented in Ref. 7, and the analog computer simulation 
developed in Ref. 8, could serve as a starting point for 
conversion to hybrid computation. 

The Naval Postgraduate School's hybrid computation 
facilities include the XDS 9300 digital computer interfaced 
with a COMCOR 5000 analog computer and an Adage Graphics 
Terminal. Perhaps the most formidable problem presented by 


a hybrid conversion would be the core storage requirement 
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in the digital computer. The present Storage capabilities 
could however be augmented by disc and/or tapes. 

The parameter scaling problems referred to in Ref. 8 
could be readily handled in a hybrid configuration. Other 
advantages to hybrid computation of the Loads and motion 
program include near real time solutions and the versatility 
of being able to change parameters during the run. The 
visual representation of the output parameters afforded 
by the Adage Graphic Terminal would be a powerful engineering 
tool. Hybrid computation would be particularly useful in 
vertical plane studies such as heave acceleration analysis 


ead prediction, and for control systems studies. 
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APPENDIX A 


An Example of Programming Changes to Improve 


Efficiency 


Reference 9 presented the following changes to the SES 


Loads and Motions Program to improve the efficiency of the 


program. The affected subroutines are RHS and INCON and in 


addition the subroutine DMINV is deleted. 


The following portion of subroutine INCON was deleted: 


Ze 


zee 


213 


214 


ZA 


DO 211 I = 1,6 

DO 211 N = 1,6 

A(I,N) = 0.0 

DO 213 N = 1,3 

A(N,N) = AM 

ACH ;4) = AIXX 

A(5,5) = AIYY 

Ge, 6) =" an ae 

Ce ee ey 

A(6,4) = -AIXZ 

AIMAX = AMAXI(AM,AIXX,AIYY,AIZZ,ABS(AIXZ) ) 
DO 214 I = 1,6 

pO 214 J = 1,6 

A(I,J) = A(I,J)/AIMAX 
CALL DMINV(A,G,D) 
POw215 1 = 196 

DO 215 J = 1,6 

AGE = A eee nx 
IF (D.NE.0.0) GO TO 10 
WRITE (6,216) 

STOP 


=) 


ae 





Inserted in place of the above programming was the 
following: 


AMASSI = 1.0/AM 
D = 1.0/(AIXX*AIZZ-AIXZ*AIXZ) 


DIXX = ALXX*D 
DIXZ = AIXZ*D 
DIZZ = AIZZ*D 
ATYYI = 1.0/AIYY 
GO TO 10 


The INCON parameters were linked to subroutine RHS by 


the following COMMON statement 
COMMON/ATRIX/AMASSI ,AIYYI,DIXX ,DIXZ ,DIZZ 


In subroutine RHS the six element matrix GF(1) was 
deleted and the following identifiers were substituted for 
the summation of forces: SUMK, SUMY, SUMZ, SUMK, SUMM, 
SUMN. In addition the following deletion was made. 

BO =a. 56 

VALUE(I) = 0.0 

pO 1J = 1,6 

VALUE(I) = VALUE(I)+A(I,J) *GF(J) 
1 CONTINUE 


cubstituted for the above DO LOOPs was the following 


VALUE(1) = SUMX*AMASSI 
VALUE(2) = SUMY*AMASSI-R*U 
VALUE (3) = SUMZ*AMASSI 
VALUE(4) = SUMK*DIZ2Z+SUMN*DIXZ 
VALUE(5) = SUMM*AIYYI 

VALUE(6) = SUMN*DIXX+SUMK*DIXZ 


The indicated changes deleted nine DO LOOPs and one 
subroutine resulting in a much more efficient program both 


from a time and a core requirement standpoint. 
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APPENDIX B 


Equivalent Expressions for Yaw Acceleration 


Reference 9 included the classic Euler equations of 
motion in six degrees of freedom and listed the assumptions 
under which certain terms of these equations could be 
disregarded as being insignificant. The question was 
raised as to whether the equation describing Yaw Accelera- 
tion (R) in classic Euler terms was equivalent to the one 
Gefined by subroutines DMINV and RHs. 

This appendix shows that, under certain assumptions, the 
two expresssions are equivalent. 

The classic Euler equations as presented in [Ref. 9! are 


as follows. 





_ F 
U = = ~ QW + RV 1) 
. +#F 
V = = Shae (2) 
P F, 
Weer FV + OU (3) 
P = Bk ZZ 4 Fn XZ o Q loz 
‘ : p 2 2 2 
em | ZG > Fs 262 HPD XZ ee es ys 
PC. eT) e 
‘ XX V¥ Viges 2 rs = CA _-TI ca XZ ) 7] 
ii afi yy ab 
Cigh Vig si 
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Qe --+—;~ + - (5) 
tyy yy tyy 
F eC) P= OR) I 

R= 7+ XX YY 4 7 XZ (6) 
ZZ ZZ ZS 


The fundamental equations of motion for six degrees of 


freedom developed in Ref. 10 are as follows: 


FL. = m[U + QW- RV -X,(Q°+R°) +¥,(PQ-R) +Z,(PR+Q)] (7) 
F, = m[V + RU- PW~Y¥,(R°+P*) + Z,(QR-P) + X_(QP+R)’] (8) 
FOS m(W + PV -QU~- Zi (PE+0*) +X, (RP-Q) + ¥,(RQtP) J (9) 
2S TP t+ (1 o-F ly) OR + m[Y,(W+PV-QU) - 2q(V+RU-PW) 
. Z eee, 
+ Xq¥_(PR-Q) - XqZ,(PQtR) + Y¥ QZ, (R -Q°~)] Olan, 
ae Ty y@t (1-12) RP + m[Z,(U+QW-RV) - XQ (W+PV-QU) 
e e 2 2 
+ ¥Q2q(QP-R) - YoX,(QRtP) + X,Zq(P aRo..] aba 
a Tight (1 y-L,,)P@ + m[X_(V+RU-PW) - Y_(U+QW-RV) 


‘ ‘ . _9 
+ ZqXqCRQ-P) - Z2,Y,(RP+Q) + ¥,X,(@ -P™)] (¥2) 


where Xa s Yq: and Za are distances from the center of gravity 


Pometemonrtela along the x, ¥, and 2 axis respectively. 
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Because of symmetry of the craft, the significant product 
terms of inertia are the XaZq Veums. — Under these conditions. 


equations (7)-(12) reduce to the following form. 


F_ = m(U + QW- RV] (13) 
Fy = aly ER = (14) 
Fo m(W+ PV - QU] (15) 
Fy = Thy? + (Iz, 71) OR - 7 (PQtR) (ilo 
Sea ey Gan ear) ees ae (17) 
m yy XX ZZ ez 
Py = I,,R+ (Iyy-T,,.)PQ + 1, (RQ-P) (18) 
where 
Lie % mXqZq 


Solving equations (13)-(18) for the desired derivatives 


yields: 
: ce 
a CN RY (19) 
} F 
v= —~-— RU + PW (20) 
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; F 


w= 2 - py + QU (21) 
EF Il _-I aL 

° kK ° 

P= > (422) ar " — (PQ+R) (22) 
XX XX XX 

: Fin Msc See Le, 232 

a ee et ae — (PO R™) (23) 
YY yY YY 

: aa i ase ee : 

R= pp - (4) PQ -  (RO-P) (24) 
ZZ ZZ ZZ 


Equations (19)-(21) are identical to equations (1)-(3). 
Upon substituting equation (24) into (22) and rearranging 
equation (23) it can be shown that equations (22)-(24) are 
equivalent to equations (4)-(6). 

If it is assumed that P, Q, and R are each much smaller 
than one, then product terms involving them may be ignored 
and equations (1)-(6) reduce to the following equations as 


stated in Ref. 9. 


: ie 
U = oy (25) 
P By 
Ve= 7 - RU (26) 
‘ F 
W = a (275) 
ijeeel Pe 
p= _k 22 = + elas : (28) 
lie 8 - I Te eee - I 
XX ZZ XZ XX Za XZ 
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F 


Q= (29) 
VV 

ae. OL ; 

RS a (30) 
ZZ ZZ 


The equations of motion as decoded from subroutines 


DMINV and RHS are as follows. 


* FY 
ye ee 
V > a oRU (32) 
; F 
: ee Ff 
P = k ZZ, ms NeeezZ (34) 
d DRT - toe = 
XX ZZ XZ XX ZZ XZ 
: is 
VY 
: lea ae geese 
R = xx 5 + es (36) 
tf -If I of -f{ | 
An aG XZ ma oe XZ 


Equations (25)-(29) are identical to equations (31)-(35). 


Substituting equation (28) into equation (30) yields: 
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R= p+ (22 dy XR Ct) ke 
ae Rug Boe XX ZZ “XZ ae 
F oi FI ¢ 
Sa se 4 _n XZ 
O° 2 2 
oe Lexie ~ tx ee ZZ -lo. XZ 
2 
FI 
ee ie. a 14 _XZ 
} no 2 o 
I I - I ZZ 
XX ZZ XZ xX EZ ZZ RT, 
2 , 2 
Jel 
ee ae 7 >, AVA VA Dy ud ZZ 7, 
7 at Bat 2 Z) ) 
T exiae ~ I xz Log ees x ) 
: eee ee 
R = + (37) 
Piece = [ Lowa ae | 
XX ZZ XZ XX ZZ XZ 


Equation (37) is identical to equation (36) therefore, based 
on the given assumptions, the classical Euler equations 

given in Ref. 9 and as developed from Ref. 10 are mathe- 
matically equivalent to those defined in subroutines DMINV 
and RHS. However it is felt that the validity of disregarding 
product terms should be more thoroughly investigated. It 

is recommended that equations of the form (1)-(6) be 
incorporated in future digital simulation studies to determine 
the degree of contribution of the cross product terms. on the 


resultant output variables. 
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